A lattice estimate of the 5'^,^^ coupling 
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We present the results of the first direct determination of the g^^,^^ coupling using lattice QCD. From our 
simulations in the quenched approximation, we obtain g^,^^ = 18.8 ± 2.3^2. J and g = 0.67 ± O.OStlj-jJg. It is in 
agreement with a recent experimental result from CLEO. 



Introduction 

The strong coupling constant g„,u^ of heavy- 
light mesons^ and a soft pion is one of the essen- 
tial parameters entering the lagrangian which en- 
globes both the chiral and the heavy quark sym- 
metry 1^ . Its precise value is useful in evaluating 
the nearest pole contribution to the shape of the 
B — > TT and D — > tt semileptonic decay form fac- 
tor, since the residua are directly proportional to 
9b'b-. and 9^,^^, respectively. As g^.^^ is not 
measurable experimentally due to phase space, its 
theoretical prediction is useful to restrict the un- 
certainties on quantities like the \Vub\ CKM ma- 
trix element. 

The coupling g„,„^ is related to the coupling 
constant 'g appearing in the chiral theory for 
heavy mesons, 

g„.B^ = 9- (1) 

The coupling g is a constant up to small 
O (l/mH) corrections in the heavy mass. A large 
number of predictions for g, using several meth- 
ods, can be found in the literature (see |^ for 
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a list of results). A surprising feature in these 
determinations is the discrepancy between QCD 
sum rules which find typically low values, g ^ 0.3, 
and quark models which naturally accommodate 
larger values^. Recently, the CLEO collaboration 
measured the coupling g^,,,-,^ and deduced, from 
eq.(|l|), a value for g (see ref Q): 

g^^^^ = 17.9 ± 0.3 ± 1.9, i.e. 

g = 0.59 ±0.07. (2) 

An exploratory lattice computation of 5 has been 
performed in ref. Q . That computation has been 
made in the static limit of HQET. To be able to 
directly confront the lattice to the experimental 
result, we decided to make the lattice QCD sim- 
ulation by using fully relativistic quarks in the 
region of the charm quark mass. Notice that, con- 
trary to the 6-quark, the charmed quarks are di- 
rectly accessible from the currently available lat- 
tices. In other words, to confront to the experi- 
ment, no extrapolation in the heavy quark mass 
is needed. 

1. Theoretical basis 

To determine 5^.^^^, we compute the matrix 
element of the axial vector current = 97/^ 75 9: 

^The Adler-Weisberger sum rule sets the bound 3 < 1. 
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between the vector [H*) and pseudoscalar {H) 
heavy-light mesons. This matrix element is 
parametrized by the three form factors Aq, 1,2(9^), 
where q is the momentum transfer between H* 
and H . 

The coupling g„,„^ is related to the form fac- 
tors through the following expression (see ref.[5] 
for details): 
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The first term (cx ^i(O)) is dominant since the 
heavy masses mu* and vtlh are almost degener- 
ate and, in the kinematical region where ~ 0, 
Ai and A2 are of the same order of magnitude. 
Our simulation has confirmed that the subdom- 
inant term (oc ^2(0)) is indeed small (< 5%). 
In studying Ai(g^), the point ^ is reached 
already at g = 0. On the lattice, it is known 
that the signals are better when no external mo- 
menta are given to the interacting hadrons. All 
these conditions make the determination of g,-,,j-,^ 
rather clean. 

2. Lattice simulation 

The main results are obtained from a simula- 
tion on a 24:^ x 64 lattice, at P — 6.2, using 100 
configurations produced in the quenched approx- 
imation with the non-perturbative 0{a) improve- 
ment for Wilson fermions. 

We have used six values of the light quark 
masses, corresponding to a pseudoscalar meson 
mass e [0.59, 0.83] GeV, and three values of the 
heavy quark masses, corresponding to rnn € 
[1.83,2.32] GeV, where, Tuh = (StriH- +mH)/4., 
is the spin averaged mass of the heavy meson. 
Note that this interval includes the value mu = 
1.97 GeV, so that no heavy mass extrapolation is 
required to reach the D-meson. 

We inserted to the current A^ the momen- 
tum q = {(0,0,0); (1,0,0)} (in 27r/L units), and 
we have checked that the case q ^ corre- 
sponds indeed to small values of q^ (we have 
g2 < 0.01 GeV^). 

The matrix elements that we compute on the 
lattice are illustrated in fig. ^ We extract the 




Figure 1. Three-point Green function computed on 
the lattice: (D\A''\D*). 
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Figure 2. Illustration of the chiral extrapolation of 
Sh-ht, ^'-"^ ^ given value of the heavy quark Wilson 
hopping parameter kq. The lower subscript "p" la- 
bels the light pseudoscalar meson. 

form factors Ai and A2 by using the standard 
analysis of the two-point and three-point Green 
functions on the lattice. 

In order to determine g^,^^ , we first have to ex- 
trapolate our simulated light quark masses down 
to the w/d-quark (corresponding to the physi- 
cal pion) and then mterpolate the heavy meson 
masses to the physical D meson. 

The chiral extrapolation has been performed 
according to several forms (linear, quadratic, 
"chiral log"; see fig. ||). In our region of light 
quark masses, we are not able to isolate non- 
linearities. Therefore, we shall use non-linear fits 
to estimate the systematic uncertainties. 

The heavy mass interpolation function is 
guided by the heavy quark symmetry, g = a + 
b/rfiH- The value of the parameter a corresponds 
to the value of g in the static limit, which we note 
doc- 

As a comment, we remark that the computa- 
tion of the form factors Ai^2, and therefore g^,^^ 
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ruH [GeV ] 


1.83(9) 2.08(10) 2.32(11) 


9 


17.7 ± 2.2 20.1 ± 2.7 22.6 ±3.3 
0.669(72) 0.668(79) 0.673(88) 



Table 1 

For each heavy-light meson mass rfiH, we show the 
values of g^,^^ and of the corresponding g. 
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Figure 3. Interpolation of g to the D-meson (hlled 
circle). 

(eq.(^)), do not depend on the improvement of 
the bare axial current (see ref. Q). 

3. Results 

We present in table 0, the values of g^'H^ 
g for each heavy-light meson mass WiH ■ All these 
quantities are already linearly extrapolated to the 
w/d-quark mass. 

The heavy quark interpolation of g to the D- 
meson mass is shown in fig. |^. Our results in- 
dicate that the slope in l/vriH for the coupling 
^ is small. Assuming that the linear dependence 
in l/rriH holds all the way to the static limit, 
l/vriH — * 0, we obtain that goo is not more than 
10 % larger than g at the level of the charm quark 
mass. 

Our final results, at /3 = 6.2, are: 



9n'D^ 



18.8±2.3l2.J, i.e 
0.67 ± 0.08l|^;;^g . 



(4) 



4. Systematic uncertainties 



In order to study the 0{a) effects, we have per- 
formed in addition to the simulation at (3 — 6.2, 



another one at /3 = 6.0. The difference between 
the two results is smaller than one standard devi- 
ation. With two values of the lattice spacing we 
cannot attempt a continuum limit extrapolation, 
but, as it is observed for other similar quantities, 
we can hope that the results at /? = 6.2 are close 
to those in the continuum limit. Further compu- 
tations at larger f3 would certainly improve the 
precision in the determination of g. 

Finite volume effects were studied by perform- 
ing simulations a,t f3 — 6.0 with two different vol- 
umes, 16'^ X 64 and 24^ x 64. Within our statistics, 
we do not see any evidence for the presence of fi- 
nite lattice volume effects. From the comparison 
between these simulations, we add a conservative 
6 % contribution to the overall systematic uncer- 
tainty. 

The chiral extrapolation is the major source of 
uncertainty in our results. We estimate this effect 
by comparing the linear fit, to the quadratic and 
"chiral log" ones (c./. Section H). 

Conclusion 

We have performed the first lattice measure- 
ment of the g^ytjj^ coupling. Our results (^ are in 
good agreement with the experimental value (^. 
The discrepancies with other theoretical predic- 
tions, like the smaller values of g given by the 
QCD sum rules, still need an explanation. Fur- 
ther improvements of our results would include, 
studying the continuum limit and attempting an 
unquenched calculation of ^. 

Work supported by the European Community's 
Human potential programme under HPRN-CT- 
2000-00145 Hadrons/LatticeQCD. 

REFERENCES 

1. R. Casalbuoni et ai, Phys. Rept. 281 (1997) 
145. 

2. D. Becirevic and A. Le Yaouanc, JHEP 9903 
(1999) 021. 

3. S. Ahmed et al. [CLEO Collaboration], Phys. 
Rev. Lett. 87 (2001) 251801. 

4. G. M. dc Divitiis et al. [UKQCD Collabora- 
tion], JHEP 9810 (1998) 010. 

5. A. Abada et at, [ ]hep-ph/0206237t ] . 



